Primordial Magnetic Fields from Dark Energy 
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Evidences indicate that the dark energy constitutes about two thirds of the critical density of the 
universe. If the dark energy is an evolving pseudo scalar field that couples to electromagnetism, a 
cosmic magnetic seed field can be produced via spinoidal instability during the formation of large- 
scale structures. 
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Recent astrophysical and cosmological observations 
such as dynamical mass, Type la supernovae (SNe), 
gravitational lensing, and cosmic microwave background 
(CMB) anisotropics, concordantly prevail a spatially flat 
universe containing a mixture of matter and a dominant 
smooth component with effective negative pressure [jjj. 
The simplest possibility for this component is a cosmo- 
logical constant. A dynamical variation calls for the exis- 
tence of dark energy whose equation of state approaches 
that of the cosmological constant at recent epochs. Many 
possibilities have been proposed to explain for the dark 
energy. Most of the dark energy models involve the dy- 
namical evolution of classical scalar fields or quintessence 
(Q). For a Q-model, the dynamics is governed by the 
scalar field potential such that the vacuum energy be- 
comes dynamically important only at recent epochs and 
gives rise to an effective cosmological constant today. So 
far, many different kinds of scalar field potentials have 
been proposed. They include pseudo Nambu-Goldstone 
boson (PNGB), inverse power law, exponential, track- 
ing oscillating, and others ||. Upcoming observations 
will measure the equation of state so as to discriminate 
between these models and distinguish them from the cos- 
mological constant || . 

Since the scalar potential V(4>) of the Q field is scarcely 
known, it is convenient to discuss the evolution of (f> 
through the equation of state, p^, = w^p^. Physically, 
1 > Wcf, > — 1, where the latter equality holds for a pure 
vacuum state. Lately some progress has been made in 
constraining the behavior of Q fields from observational 
data. A combined large scale structure (LSS), SNe, and 
CMB analysis has set an upper limit on Q-models with 
a constant w^, < — 0.7 (J], and a more recent analysis 
of CMB observations gives = -0.82±£^ §. Fur- 
thermore, the SNe data and measurements of the posi- 
tion of the acoustic peaks in the CMB anisotropy spec- 
trum have been used to put a constraint on the present 
< —0.96 Q.The apparent brightness of the farthest 
SN observed to date, SN1997ff at redshift z ~ 1.7, is 
consistent with that expected in the decelerating phase 
of the flat ACDM model with Q A ~ 0.7 @, inferring 
= —1 for z < 1.7. In addition, several attempts have 



been made to test different Q-models [B. Nevertheless, 
it is primitive to differentiate between the variations, and 
the reconstruction of V((j>) would require next-generation 
observations. 

Although the physical state of the dark energy can be 
probed through its gravitational effects on the cosmolog- 
ical evolution, it is important in fundamental physics to 
understand whether the quintessence is a nearly massless, 
slowing rolling scalar field. It has been pointed out by 
Carroll that the existence of an approximate global sym- 
metry would allow a coupling of the Q field, <j>, to the 
pseudoscalar F^F^ of electromagnetism, which pro- 
vides a potential observable in polarization studies of 
distant radio sources |J. As long as an ultralight 4> ne ld 
couples to photon where for a slow-roll condition, the 
mass is comparable to Hq, it is conceivable to have 
very long-wavelength electromagnetic fields generated via 
spinodal instabilities from the dynamics of <f> as a possible 
source of seed magnetic fields for the galactic dynamo. 

As we know, the issue of the origin of the observed clus- 
ter and galactic magnetic fields of about a few pG |l(]] 
remains a puzzle jy]. These magnetic fields could have 
been resulted from the amplification of a seed field of 
B S eed ~ 10 _23 G on a comoving scale larger than Mpc via 
the so-called galactic dynamo effect. A number of scenar- 
ios have been proposed for generating seed fields in the 
early universe, mainly relying on non-equilibrium con- 
ditions such as inflation, the electroweak and the QCD 
phase transitions. After the phase transitions, the uni- 
verse became a highly conducting plasma so that the 
magnetic flux which existed is frozen in, and the ratio of 
the magnetic energy density and the thermal background, 
Pb/Pj, remains constant thereafter. The required B see d 
translates into ps — 10~ 34 p 7 . However, it turns out that 
the generated fields in these models are too small, except 
in somewhat contrived cases, to be of cosmological inter- 
est. In an equilibrium condition, a large damping term 
induced by the high plasma conductivity suppresses sig- 
nificantly any electromagnetic field fluctuations jl2],[L3| . 

In this Letter, we will investigate the implication of 
the electromagnetic coupling of the evolving </> field to 
the origin of the primordial magnetic field (PMF). Then, 
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we will provide a mechanism in which the PMF can be 
generated via the 0-photon resonant conversion during 
the LSS formation. 

Here we consider the 0-photon coupling, 



L^-y = -0 < 



F a pF t 



(1) 



where — d^Ay ~ d v A^, / is a mass scale, and c is a 
coupling constant which we treat as a free parameter. For 
the present consideration, we pick / equal to the reduced 
Planck mass M p = (87rG) -1 / 2 . We are thus led to study 
the cosmological evolution of the (/)-photon system in a 
flat universe. The effective action is 
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where the signature is ( — h ++) and Sm is the classical 
action for matter. We assume that the universe today 
has matter il M = 0.3 and quintessence 175 = 0.7, and 
define an f^-weighted average jl] 
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where 770 and rji s are respectively the conformal time 
today and at the last scattering, defined by rj = 
Hq J cfta _1 (i) with the scale factor a and the Hubble con- 
stant Hq. Assuming a spatially homogeneous <p field, the 
equation of motion is given by 
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where 9 — <j>/M p . The last term of Eq. (|J) is the back 
reaction from the produced magnetic fields. Later we will 
show that this term is too small to have any effect on 
the evolution of the 4> field when we consider the photon 
production. For now we omit this term and thus the 
evolution of the cosmic background is governed by 
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where h = H/H , = p^j '(M p H ) 2 , and we have used 
<j> 2 = (1 + w^Pa and V((j)) = (1 — w < f,)p t j > /2. We have nu- 
merically solved the background equations by proposing 
a simple square-wave form for w,p as shown in Fig. |l|. In 
order to satisfy the above-mentioned observational con- 
straints on Wa,, we have chosen = — 1 for z < 2, and 
a width of the square-wave such that (uu) ~ —0.7. Note 



that at the present time fl^ = Pc /,/{3h 2 ) = 0.7 and H t = 
0.95. We have also plotted d9/drj — a^J (1 + wa)pa- 

In the original PNGB models for dark energy , the 
well-known PNGB cosine potential was used. This re- 
sults in that > at the present time, which is dis- 
favored by the observations. In fact, the square-wave 
equation of state is anticipated and quite general in the 
class of Q-models using PNGB fields incorporated with 
quantum corrections to the cosine potential |l7| . In 
the PNGB models, the nonperturbatively large quan- 
tum fluctuations driven by spinoidal instabilities strongly 
dissipate the oscillation of the scalar field expectation 
value, leading to the field oscillation to become over- 
damped, and finally the field expectation value relaxes 
to the minimum of the effective action. However, we 
would like to emphasize here that the spinodal instabil- 
ities that we will discuss later to drive the amplitude 
fluctuations of long-wavelength magnetic fields to grow 
is a general phenomenon which is not restricted to such 
a square-wave equation of state. In fact, in the present 
consideration, the spinoidal instabilities responsible for 
the generation of the PMF occur during the formation 
of LSS for z ~ 10 — 2 as long as during which the equa- 
tion of state of the scalar field deviates from = — 1. 
After z ~ 2, the quintessential potential dominates and 
behaves like a cosmological constant with w ~ — 1 that 
ceases the photon production. Several quintessence mod- 
els have been proposed to generate the equation of state 
of the scalar field that has the above-mentioned features. 
In Ref. @, a potential of the form V{4>) = Vpi^e'^, 



where V^W ^ s a polynomial of (f> and A is a parameter, 
has been invented to introduce a local minimum in the 
exponential such that the field gets trapped into it. After 
the field gets trapped it starts behaving like a cosmologi- 
cal constant and the universe eventually enters an era of 
accelerated expansion. With suitable parameters, they 
showed that the equation of state stays near the value 
of about 0.3 for high redshifts and abruptly decreases to 

— 1 for z < 3. This desired accelerated expansion can 
be also provided by a different V p (<p) arising from the in- 
teraction between branes in extra-dimension models [ p"4"| , 
or having an oscillating term in the tracking oscillating 
model by Dodelson et al. in Ref. Q. In fc-essence mod- 
els with a modified kinetic term for the scalar field, the 
equation of state suddenly drops to the value of about 

— 1 at z ~ 2 Jl5]| . In summary, we have reviewed some 
quintessence models which have the desired evolution of 
the scalar field. With or without further tuning of the 
model parameters, they can be used to provide a cosmo- 
logical background for the generation of the PMF during 
the LSS formation with the photon production ending at 
z ~ 3. However, for our purpose to illustrate how the 
PMF can be generated from the quintessence dynamics, 
we will use this simple, generic model in Fig. |l| that sat- 
isfies all existing observational constraints to carry out 
the calculation. 

The most stringent limit on the (/)-photon coupling 
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comes from the cooling via the Primakoff conversion of 
horizontal branch (HB) stars in globular clusters [fL8| , 
c/f < 1.5 x Kr n GeV-\ which gives c < 3.7 x 10 7 . If 4> 
carries a mass ~ Hq , it would decay into two photons 
with a width, T = c 2 m^/ (Airf 2 ). Hence the lifetime of 
(j> is much longer than the age of the universe. There is 
a very nice limit on c coming from the rotation of the 
plane of polarization of light from distant radio sources, 
c < 3 x 1O _2 /|A0|, where A9 is the change in 9 between 
z = 0.425 and today ||. Since we have taken = — 1 
(i.e., A9 — 0) for z < 2, this limit does not apply to 
our case at all. From the theoretical point of view, it 
has been argued that because of radiative corrections c 
must be infinitesimally small to keep the <f> field light [jl9| . 
However, in this study we would like to find out the value 
of the 0-photon coupling needed in order to produce the 
large PMF over large correlation length scales if we as- 
sume that the generation of the PMF is solely due to the 
quintessence dynamics during the LSS formation. 

From the action (|^), the comoving magnetic field in 
the comoving coordinates (r, x) with r = rj/Ho satisfies 
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where V = d/dx, a is the plasma conductivity, and 
v = dx/dr is the peculiar plasma velocity. After hy- 
drogen recombination, the residual ionization keeps the 
conductivity high, a/H ~ 10 22 (T/eV)~ 3/2 > 1. As a 
result, the a term under the assumption that v = 
damps out any growth of the B field on scales above 
~ A.U. flj2|| . But this assumption may not hold when 
the universe has entered the non-linear regime. It could 
be understood from the recent study in Ref. J2(J where 
the battery mechanism has been proposed as a source 
for generating a small initial magnetic field during LSS 
formation. 

From the theory of small fluctuations, v will grow with 
density perturbations. Velocity flows grow rapidly af- 
ter extreme nonlinearities develop in the cosmic fluid. 
Therefore, the authors in Ref. |2Q] performed numerical 
simulations by adding the battery equation to the hydro- 
dynamic code for structure formation, and showed that 
magnetic fields are built up in regions about to collapse 
into galaxies. The non-linear turbulence effect plays an 
important role in generating and subsequently amplifying 
the created PMF which overcomes the damping mecha- 
nism due to the high plasma conductivity. The numerical 
simulations showed that the rms value for V- v ~ V x v is 
about 10~ 16 s~ 1 on comoving scale of Mpc at z ~ 3 [ po[ . 
This is just the inverse of the time scale Tlss f° r LSS 
formation. When a Mpc-scale magnetic field with mag- 
nitude B is simultaneously created, and as long as the 
growth rate dB/dr ~ B/tlss is comparable to the twist- 
ing term | V x (v x B)| ~ BV ■ v, the a term would be sig- 
nificantly reduced and the high plasma conductivity is no 



longer a hindrance to the growth of large-scale magnetic 
fields. Here we will show that the temporal variation of 4> 
can also be a generating source for large-scale magnetic 
fields. However, we will not pursue numerical simula- 
tions similar to Ref. p0| , combining their battery source 
with the quintessence dynamics. As a first step, we will 
simply omit the a term and solve for the photon equa- 
tion self-consistently, while following closely their results 
and estimating the effect of the magneto- hydrodynamical 
dissipation. Of course, it would be very interesting to 
carry out direct numerical simulations by coupling the 
</>-photon system in Eq. (Q) to the hydrodynamic code to 
make a map of the created magnetic fields. 

Now we write B = V x At and decompose the trans- 
verse field At(t, x) into Fourier modes, 
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where 6± k are destruction operators, and e± k are circular 
polarization unit vectors. Then, defining q = k/Hg, the 
mode equations are 



drj 



with initial conditions at early epoch given by 

v ±q = i, = -*<?■ 
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Hence, the comoving energy density of the magnetic field 
is given by pb = (B 2 )/8tt = j dq{dps / dq) with 
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where the coth term is the Bose-Einstein enhancement 
factor due to the presence of the CMB with current tem- 
perature To and energy density p 1 — 7t 2 Tq/15. The 
nonlinear growth in V+ q can be understood by treating 
the background field solution d9jdr\ > in Eq. (B) as a 
constant. For small fluctuations of the field, the long- 
wavelength modes of V+ q , where q lies within the unsta- 
ble band, q < Ac(d9/dri) ) in fact "see" the inverted har- 
monic oscillators that cause the amplitude fluctuations 
to grow exponentially, while the evolution of V- q modes 
is purely oscillating without any growth. The growth 
of these unstable modes is driven by spinodal instabili- 
ties pl| . Now we can write the approximate solution for 
the growing modes to Eq. (^) as V+ q cx e av where the 
exponent a — Acq(d9 / drj) — q 2 provided that q is in 
the unstable band. The maximum value of a is a max = 
2c(d9/drj) at q = 2c{d6/drf). Since the maximum growth 
for the amplitude fluctuations occurs at q — 2c(d9/drj), 
in order to produce the large-scale correlated magnetic 
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fields on 10 Mpc scales with the fastest growth, one can 
estimate the order-of-magnitude of the coupling c, which 
is about 10 2 , using the fact that dO/d-q ~ 0.8 at redshift 
z = 10 after which the universe was in the nonlinear 
regime and based on our previous argument that the high 
plasma conductivity effect which could damp out the pro- 
duced magnetic fields can be consistently ignored. The 
self-consistency condition to reduce significantly this high 
conductivity effect can be justified in the sense that for 
the exponentially growing modes that we are interested 
in, dB/dr ~ a max H B ~ (lO^s" 1 ^ with the cou- 
pling c ~ 100 and d9 jdr\ ~ 1, which is of the same order 
of the twisting term |V x (v x B)| ~ BV ■ v during the 
LSS formation obtained from the numerical simulations 

As you will see in our numerical analysis, the spinodal 
instability provides a robust mechanism to generate the 
long- wavelength fluctuations where the created magnetic 
fields correspond to coherent collective behavior in which 
the fields correlate themselves over a large distance of or- 
der of 10 Mpc. However, since d9/drj is actually time 
dependent, the modes will grow only for a period of time 
before they move out of the unstable band. In the end, 
the growth of the long-wavelength modes will be shut 
off completely when the bandwidth of the unstable band 
vanishes as d6 jdr\ reaches zero at redshift z ~ 4. We have 
numerically solved the mode equations (^|) using c = 130 
and the background solution as shown in Fig. U and plot- 
ted the ratio (dps / dq) / p 7 in Fig. |^. Although photons 
are being produced as the scalar field starts rolling at 
z ~ 60, we have counted the photons produced only af- 
ter z = 10 when the universe has presumably entered the 
non-linear regime. The result shows that a sufficiently 
large seed magnetic field of 10 Mpc scale has been pro- 
duced before z ~ 4. Moreover, we notice that when 
c ~ 130 the spinoidal instability gives rise to magnetic 
seed fields of the right magnitude and length scale. If we 
decrease the value of c, the peaks of the curves in Fig. ^ 
will shift to the lower left-hand corner. 

The magnetic energy density can be approximately ob- 
tained from Fig. | and is about fi s ~ 10~ 33 < Q.^. 
Therefore, we anticipate that the back reaction effect 
from the produced magnetic fields on the scalar field evo- 
lution equation is negligible. This can be shown by eval- 
uating the last term of Eq. with 
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which is found to be extremely small compared to the 
other terms in Eq. (Q) when the photons are being pro- 
fusely produced. 

Although the value of c = 130 is well below the HB 
limit, undoubtedly it is much larger than the theoreti- 
cal expectation. However, as suggested in Ref. ||, an 
unsuppressed </>-photon coupling may arise in higher di- 
mensional theories. One possible way to reduce the value 



of the coupling and at the same time to produce the suffi- 
ciently large PMF on large correlation length scales dur- 
ing the LSS formation is to combine the mechanism we 
proposed here with the battery source in Ref. fl2Q| . Also, 
note that the growth rate is controlled by the exponent 
a ~ 2c{d9/drj). As such, it may be possible to reduce the 
value of c by having a large dO/drj. In the case of scalar 
quintessence, we have seen in Fig. [j] that d9/dr] < 0.8 for 
z < 10. We cannot further increase d9/drj since has 
already reached the maximum value. Perhaps, in the k- 
essence models, in which the kinetic term is modified to 
K((j))(p 2 /2 and we have <j) 2 — (1 + w^)p^j K ((f)), we may 
make d9/dr] much larger than one by tuning the prefactor 
K(4>). This possibility is under investigation. 

In conclusion, we have made an effort to link the dark- 
energy problem to a solution to another important prob- 
lem in cosmology, namely, the generation of primordial 
magnetic fields. So far, we can only probe the equation of 
state of the dark energy through its cosmological effects. 
Undoubtedly, it is extremely important to have any clue 
about the nature of the dark energy. Here we have stud- 
ied the possibility of generating the PMF on Mpc scales 
during the LSS formation by coupling the electromag- 
netic field to the evolving scalar field that accounts for 
the dark energy dynamics. The high conductivity effect 
due to residual ionization after hydrogen recombination 
can be argued to be significantly reduced as a result of 
the existence of the cosmic flow with nonlinear, twist- 
ing peculiar velocity to avoid a hindrance to the growth 
of the magnetic fields. We have shown that the nonlin- 
ear instability that drives the rapid growth of magnetic 
fields is of spinodal instability where the long-wavelength 
modes about the Mpc scales evolve as being the inverted 
harmonic oscillators and the amplitude fluctuations be- 
gin to grow up to a time at which the scalar field velocity 
approaches zero at redshift z ~ 4. We found that when 
c ~ 130 the spinoidal instability gives rise to magnetic 
seed fields of the right magnitude and length scale. This 
strong coupling is well below the HB limit, but much 
larger than the theoretical expectation. 

Here we have used a standard scalar quintessence and 
a simple equation of state as a working model. It is inter- 
esting to extend the idea to non-standard scalar models 
such as A;-essence to see if the required 0-photon cou- 
pling can be reduced to a smaller value with the specific 
equation of state. On the other hand, new laboratory 
searches for the coupling of photon to pseudo-scalar and 
more ingenious astrophysical or cosmological limits on 
the coupling would be needed. In addition, future CMB 
and SNe experiments would put a strong constraint on 
the quintessential equation of state at high-redshifts. In 
particular, combining the battery mechanism in Ref. PQ] 
with the 0-photon coupling which provides an alternative 
robust mechanism to generate the PMF in the hydrody- 
namical simulations would be a very interesting subject 
to tackle. 
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FIG. 1. The quantities d6/drj, W4,, Hot, and fi^, as a func- 
tion of redshift. Note that d6/drj is drawn 4 times smaller. 




FIG. 2. Ratios of the spectral magnetic energy density to 
the present CMB energy density at various redshifts. The 
present wavelength of the magnetic field is given by 2tv / (qHo) ■ 
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